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ABSTRACT 

Simultaneous measurement of line-of-sight (LOS) magnetic and velocity fields 
at the photosphere and chromosphere are presented. Fe I line at A6569 and at 
A6563 are used respectively for deriving the physical parameters at photospheric 
and chromospheric heights. The LOS magnetic field obtained through the center- 
of-gravity method show a linear relation between photospheric and chromospheric 
field for field strengths less than 700 G. But in strong field regions, the LOS 
magnetic field values derived from are much weaker than what one gets from 
the linear relationship and also from those expected from the extrapolation of the 
photospheric magnetic field. We discuss in detail the properties of magnetic field 
observed in from the point of view of observed velocity gradients. The bisector 
analysis of Stokes / profiles show larger velocity gradients in those places 
where strong photospheric magnetic fields are observed. These observations may 
support the view that the stronger fields diverge faster with height compared to 
weaker fields. 



Subject headings: Sun: magnetic fields - Sun: photosphere - Sun: chromosphere 
- sunspots 
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INTRODUCTION 



It is inferred through various observations that the magnetic field is playing a central 
role in the solar energet i cs tha t take place in the higher layers of the atmosphere (see for 



eg. iRegnier fc Canfieldl (l2006l )). But, it is somewhat difficult to obtain a reliable vector 
magn etic field measurements at the chromosphere and corona (ISocas-Navarrd l2005l : [Judge 
20071 ). The measurements of the photospheric magnetic fields are relatively well established. 
However, simultaneous vector magnetic field observations at the photospheric and chromo- 
spheric heights will give a better handle on the understanding of the magnetic structuring 
of the solar atmosphere. More direct and reliable measurements of the magnetic field at the 
chromosphere will also serve as a better boundary condition for extrapolating the magnetic 
fields to the coronal heights. 

Measurements of the vector magnetic field in Ha is particularly important in under- 
standing the connection between photospheric, chromospheric and coronal magnetic field as 
i ts height of formatio n ranges almost from the upper photosphere to upper chromosphere 
( IVernazza et al.lll98ll ). Also, this is one of th e most widely used spectral lines for the study 
of solar chromosphere (see a recent review by iRuttenl (120071 ) ) . 



Comparison of the active region magnetic fields measured in and in the photosphere 
show a one-to-one correspondence in the weak field regions. But in the strong field regions. 



like umbra, there is a considerable deviation from the linearity (IBalasubramaniam et al. 



2004 : lHanaokall2005l ). Infact the LOS magnetic field measured in Ha weakens much faster 
for the corresponding strong field regions of the photosphere. At this point it may not be 
appropriate to demarcate the field strength above whi ch this deviation happens b ecause 



different observations show different deviation points. iGosain &: Choudharyl (120031 ) have 
also reported the systematic weakening of the magnetic field derived from Mg I A5173/5184 
lines in comparison with that of Fe I A6301. 5/6302. 5 lines. In their case, the magnetic field 
measured by Mg I lines agree with the potential field extrapolation of the photospheric LOS 
magnetic field in weak field regime. In strong field regions there is a systematic shift towards 
lower values but still linear. Simultaneous observations of He I at A10830 and Si I lines at 



A10827.1 by lChoudhary et al.l (120021 ) suggest that the field diverges faster in the upper layers 
of the chromosphere. 

However there is no conclusive explanation for the observed weaker chromospheric 
magnetic fields fo r the corresponding strong pho tos pheric fields. Va rious possibilities have 
been discussed by iBalasubramaniarn et al. ( 2004) fclHanaokal ( 20051 ) about weaker chromo- 
spheric fields observed in umbra. iHanaokal (120051 ) have discussed the possibility of scat- 
tered light and/or peculiarity of the atmosphere (radiative transfer effects) for the decrease 
in the polarization signal which in turn cause the underestimation of the magnetic field. 
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Balasubramaniam et al.l (120041 ) have suggested that the strongest fields measured at the 
photosphere diverge spatially and more quickly than the weak fields when propagating up- 
ward. In this paper, we address these issues from the point of view of the observed velocity 
gradients. 

It is well known th at, velocities and magn etic fields are coupled to each other in the solar 
atmosphere (see for eg. iRajaguru et al.l (120061 )). Any change in the magnetic field is expected 
to alter the plasma motion and hence the observed velocity. If the field strength decreases 
with height then the velocity is expected to increase and vice versa, for the simple reason 
that the inhibition of the plasma depends on the magnetic field strength (jSpruit fc Zweibel 
19791 ). It is also well known that the magnetic field strength in the umbral region at the 
photosphere is large. If the magnetic field measured with Ha in the umbral region is weaker, 
then it shows that the magnetic field gradient is larger than expected which implies larger 
velocity gradients in the umbral region compared to the penumbral region at chromospheric 
heights. In order to observationally verify this we have analyzed the velocity and magnetic 
fields estimated from the simultaneous spectropolarimetric observations of an active region 
using Ha and Fe I spectral lines. 

The spectropolarimetric observations which were carried out at the Kodaikanal solar 
telescope using a dual beam polarimeter are discussed in section([2]). Sections (3) and (4) 
discuss the data reduction and analysis procedures respectively. A comparison of LOS mag- 
netic field at the chromosphere and photosphere is presented in section (15.11) . Velocity 
gradients obtained through bisector technique are discussed in section (15.21) . 



INSTRUMENT AND OBSERVATIONS 



Spectropolarimetric observations were carried out usin g a newly added dual beam po- 
lari meter t o the s pectrograph at Kodaikanal solar telescope ( iNagaraju et al.l (120071 ) and also 



see iBappij ( 119671 ) for details about the spectrograph and telescope setup). The wavelength 
region of the observations presented in this paper includes if„(A6563) and Fe I (A6569) lines. 
Both Hn and Fe I lines are magneti cally sens itive with effective L ande 'g' factor of 1.048 
( ICasini fc Landi Degrinnocentilll994l ) and 1.4 (IKobanov et al.ll2003l ) respectively. The mea- 
sured instrumental broadening and the linear dispersion close to this wavelength region in 
second order of diffraction are 38 ± 0.5mA and 10.15mA/ pixel respectively. Th e wavelength 
calibration was done using the telluric lines {H2O) at A6570.63 and A6561.097 (IMoore et al. 
1966h . 



An eight stage modulation scheme was used for the measurement of the general state 
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of polarization (iNagaraju et al.l 120071 ) . In this modulation scheme the measurement of the 
Stokes parameters are well balanced over the duration of eight stages of intensity measure- 
ments. The rotation of the retarders were done manually for the modulation of the input 
light. 

The spectropolarimetric data of an active region NOAA0875 presented in this paper 
were obtained on 28th April, 2006. The heliographic coordinates of the sunspot during 
observations were 11° south and 18° west. Scanning of the sunspot was done by moving the 
Sun's image in an east- west direction in steps of ~ 5". For each sht position, the modulated 
intensities were recorded by the CCD detector. The CCD is a IK x IK Photometries detector 
with the pixel size of 24/i. Eight stages of modulation took about 90 s with a typical exposure 
time of 0.5 s. 



DATA REDUCTION 



The spectral images were corrected for dark current and gain table v ariation of the pixels 

(detai ls about the flat field ing of spectropolarimetric data can be found in lSchlichenmaier fc CoUados 
fl2002h : lBeck et al.l toO^ ). 



Sankarasubramanian fc Rimm 



Model independent velocity and magnetic fleld gradients can be derived from the bisec- 
tor technique for the range of heights over wh ich the spectral line is formed, provided the full 
spect ral proflle is available without any blend ( Balasubramaniam et al.lll997 
20021 ). The apphcation of bisector technique to Ha is restricted because of the blend in its red 
wing by the Zeeman sensitive Co I (A6563.4) line which forms at the photospheric heights. 
In the following section we discuss the procedure to remove this blend and its limitations. 



3.1. Blend Removal Procedure 

Individual proflles of Ha corresponding to the orthogonally polarized beams in each 
stage of modulation were considered for removing the blend. A function which is a linear 
combination of a Gaussian and a quadratic term was fltted to the blend region of the observed 
Ha proflle. This non-linear least square flt (available in IDL) takes into account the curvature 
in the intensity proflle of Ha line along with the Co I line proflle approximated as a Gaussian. 
The Gaussian, constructed out of the fltted parameters, was removed from the observed 
proflle. Then the eight intensity measurements were combined to obt ain Stokes /, Q, U an d 



V spectral images through the demodulation procedure explained in INagaraju et al.l (120071 ). 



A typical Ha intensity proflle before and after the blend removal is shown in Fig. [T]with the 
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bisectors marked as diamond symbols. 

Even though it appears in the total intensity profile that the effect of Co I line is 
completely removed, the blend residuals still appear in Stokes Q, U, and V profiles. The 
residuals are due to the Gaussian approximation used for the Co I line profile. However, 
it is demonstrated in the section 15.21 that the blend residuals do not have any effect on the 
velocity gradients calculated from the Ha intensity profiles. 

The better way to remove the blend may be to synthesise the Stokes profiles of Co I 
using radiative transfer equations along with the atmospheric parameters obtained through 
Fe I (A6569) line. However in this paper only Stokes / profiles are considered for the velocity 
gradients estimation and a restricted spectral range of Ha about line center for the LOS 
magnetic field estimation. Hence, the simple blend removal method outlined in this paper 
is found to be sufficient. 



3.2. Correction for Polarimeter Response and Telescope Induced Cross-talks 



The polarimetric data was corrected for polarimeter response (for details see lNagaraju et al. 



(120071 )). The instrumental polariz ation introduced by the t elesco pe w as corrected by using 
the te lescope model developed by iBalasubramaniam et al.l (119851 ) and ISankarasubramanian 
(120001 ). The refra ctive index va lues for Aluminium coating used in the model are obtained 
from the catalog (IWalter Ill978l ). Since these values may be different from the actual values, 
there still remain residual cross-talks among S tokes parameters. These re s idual cross-talks 
are removed by the statistical method given in ISchlichenmaier fc CoUadod (120021 ). 



4. DATA ANALYSIS 



4.1. LOS Magnetic Field 



The LOS magnetic fields at the photospheric and chromospheric heights are derived 
using the Fe I line and Ha respectiv ely. For the derivation of the LOS magnetic fields, center 



of-gravity (COG) method was used (IRees fc Semellll979l : ICauzzi et al.lll993l : IUitenbroekll2003 
Balasubramaniam et al.ll2004l ). 



For the COG method, the LOS field strength is given by 



B 



LOS 



(A+ - A_)/2 
4.667 X 10-13 Ag^i 
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where Aq is the central wavelength of the line in A, ql is the Lande 'g' factor of the 
line, and X± are the COG wavelengths of the positive and negative circularly polarized 
components respectively. The COG components are calculated as 

where Icont is the local continuum intensity. The integration is over the spectral range 
of a given spectral line. Since Fe I line is free of blend, the full spectral range is available 
for the COG method. For the Ha line, the spectral range is restricted due to the blending 
of the Zeeman sensitive Co I line. Even though the blend is removed in Stokes / profile 
using the Gaussian fit technique explained in the section (13.11) . there still remains a residual 
in polarization profiles. To avoid introducing any artifacts in the LOS magnetic field values 
derived using COG method, a restricted spectral range about the Ha line core is used. This 
spectral region was selected by looking at the strongest Stokes V signal of Co I line. Hence 
the derived LO S magnetic field using Ha would correspond mostly to the upper chromosphere 



f lRuttenl 120071 ) 



The expected maximum underestimation of the photopsheric LOS magnetic field is 
about 12% for the kind of field strengths presented in this paper. Because of the large 
intrinsic Doppler width of Ha, the underestimation of the LO S field is not exp ected even 



for the strongest magnetic field observed at the chromosphere ( lUitenbroekl l2003l ) . However 
we checked the reliability of the COG method in deriving LOS field from Ha by applying 
it to synthetic St okes I and V profiles of Ha obtained using the radiative transfer code of 



Uitenbroekl (119981 ). From these studies it was found that the COG method overestimates the 
LOS field by about 1.5% in the field strength range to 2000 G. However, the main error 
involved in the determination of the LOS magnetic field, for this observation, is due to the 
error in estimating the shifts between the COG wavelength positions X^ and A_ which is 
expected to be less than 70 G. 



4.2. Photospheric Vector Magnetic Field 



The strength and the orientation of the chro mospheric magneti c field are very much 
dependent on its vector nature at the photosphere (jWiegelmann et al.l (120061 ) and references 
there in). Inference of the vector magnetic field at the chromosphere from the Ha observa- 
tions is difficult because. Stokes Q and U are hardly discernible. Only at few locations of 
penumbra they are above the noise level (2 x 10~^ I^ont in our observations). Hence for the 
correct interpretation of the LOS magnetic field obtained from Ha it is important to have 
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the complete information about the vector magnetic field at the photosphere. 

The photospheric vector magnetic fields are obtained by inverting the observed Stokes 
profiles of Fe I line. Milne-Eddington Line Analysis using a Numerical Inversion Engine 
(MELANIE) was used to perform the inversion. MELANIE performs non-linear least- 
square fitting of the observed Stokes profiles under Local- Thermodynamic-Equilibrium (LTE) 
condition by assuming Milne-Eddington atmosphere. Inversion code returns magnetic field 
strength, inclination angle with respect to LOS, azimuth, line strength, damping parame- 
ter, LOS velocity, source function and its gradient with optical depth, macroturbulence and 
fraction of stray light/fill factor of the non-magnetic component. 

The fit error in estimating the magnetic field ranges from 50 G for the Stokes profiles 
which are symmetric and well above the noise level to 250 G for the Stokes profiles which 
are highly asymmetric and close to noise level. The maximum fit error in the estimation of 
the field orientation is about 5° and its azimuth is 6°. Other physical parameters returned 
by the MELANIE are not used in the current work and hence are not discussed here. 



4.3. Velocity Gradients 



The measuremen t of the magnetic field in has been difficult to interpret (IBalasubramaniam et al 



2004 : lHanaokall2005l ). However it is possible to address some of these difficulties from the 
study of plasma m otions asmagnetic field and plasma motion infiuence one another in the 
solar atmosphere (IGaryl 1200 ll ). Hence the study of velocity and its gradient may help in 
better interpretation of the observed magnetic field. 

The velocities at the photosphere and chromosphere are obtained through COG method 
( lUitenbroekl l2003l ) . The COG wavelength Xcog of a fine profile / is defined as the centroid 
of its residual intensity profile: 



A 



COG 



Jihont - I)d\ ' 



(3) 



The LOS velocity with respect to the average quiet Sun reference (A^e/) is defined as 

c(AcoG " Kef) 



vlos 



X. 



ref 



(4) 



where c is the speed of the light. 



http://www.hao.ucar.edu/public/research/cic/index.html 
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Bisector technique has been apphed to Stokes / profiles to derive the velocity gradients 
both at the photosphere and chromosphere. The Stokes / profiles of Fe I and are 
separately considered for bisector analysis. Bisectors are obtained at 9 equal intensity levels 
between line core and the wing for Fe I and 14 equal intensity levels between line core and 
wing for Ha respectively. Out of 9 bisectors of Fe I line only 7 are considered, namely, 
the bisectors between second and eighth counting from the line core. Similarly for Ha the 
bisectors between second and thirteenth are considered totalling 12 bisectors. The bisectors 
very close to line core and wing are not considered because of the lower signal-to-noise ratio 
and to avoid the influence of the continuum respectively. For the Fe I line, the wavelength 
position of the second bisector which corresponds to higher atmospheric layer was subtracted 
from the seventh bisector which corresponds to lower atmospheric layer. Similarly for Ha-, 
the wavelength position of the second bisector was subtracted from the thirteenth bisector. 
The wavelength differences (AAs) thus obtained are converted into velocities - which would 
then represent the velocity difference between the lower and higher atmospheric layers - using 
the following relation, 

cAA 

AHs = (5) 

Where Aq is the rest wavelength of the spectral line under consideration. The velocity 
difference defined in Eq. (E]) represents the velocity gradient over the line formation height. 

The errors in estimating the velocity differences are mainly due to the errors involved 
in finding the wavelength shifts. The maximum error in estimating the velocity gradients is 
about 0.09 km s~^. 



4.4. Stokes V Amplitude Asymmetry 

Amplitude and area asymmetries of Stokes V profiles are caused b y the gradients in the 

veloci ty and magnetic fields (see for eg. ISanchez Almeida fc Lited ( 1l992l ) ; ISankarasubramanian fc Rimmele 
( I2OO2I )). Hence their analysis will give us some handle on the understanding of the field gra- 
dients. 



If ar and aj, represent the amplitudes of red and blue wings of Stokes V respectively 
then the amplitude asymmetry is defined as. 



5a 







Or 


Ob 


+ 





(6) 



The area asymmetry is not considered in this paper due to the difficulty in estimating 
it for the Ha in the presence of Co I line blend. 
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5. RESULTS AND DISCUSSIONS 

5.1. Comparison of Photospheric and Chromospheric LOS Magnetic Fields 

The scatter plot of LOS magnetic field derived from and Fe 1 is shown in Fig. 
[21 The plot shows that the chromospheric magnetic field is weaker in general compared 
to its photospheric counterpart. The important point to note from this figure is that the 
chromospheric fields are much weaker in the locations where the s trong photospheric fields 



are ob served. Similar kind of observations were reported earlier by iBalasubramaniam et al. 



( 20041 ) using simultaneous observations in and Fe I line at A6301.5 and i Hanaokal J2005h 



who compared the LOS field measured in Ha with the magnetograms of SOHO/MDI. These 
observations may imply that the stronger fields weaken much faster when they propagate 
upward in the solar atmosphere. 

To illustrate the quicker weakening of the stronger fields, plots of the photospheric and 
the chromospheric LOS field strengths along two different radial slices of the sunspot are 
shown on the right side of the Fig. [31 The radial slices considered for this purpose are 
marked as 1 and 2 on the SOHO/MDI intensitygram showing the sunspot analysed in this 
paper B The top panel on the right side of the Fig. [3] shows the plots of photospheric and 
chomospheric LOS magnetic field strengths along the radial cut (marked as 1 on the sunspot 
image) passing close to the central umbra. Note that, along this radial cut the photospheric 
field strength increases systematically from the penumbral region to the umbral region. 
While the chromospheric field strength increases upto about 800 G, more or less linear with 
the photospheric field, and then starts decreasing towards the umbra. The decrease in the 
chromospheric field is much larger close to central umbra. The photospheric field along the 
radial cut 2 as shown in the bottom panel on the right side of the Fig. [31 shows the behavior 
very similar to that was seen along the radial cut 1. That means the photospheric field 
strength is larger in the umbral region and decreases towards the edges of the sunspot. In 
the case of chromospheric field strength, the values do not decrease towards the umbra as 
was seen for the radial cut 1 but, they are considerably smaller compared to its photospheric 
counter part. The main difference between these two radial cuts is that the field strength 
in the umbral photosphere for radial cut 1 is larger compared to that of radial cut 2. This 
may be an indication of the faster divergence of the stronger fields. However, there are other 
possibilities which can cause this observed weaker fields in the umbral chromosphere and 



^Because of the coarser step size used for scanning the sunspot (w 5 arc sec), the raster images do not 
give a good representation of the observed region and also there was no imaging facihty available during 
these observations and therefore intensitygram from SOHO/MDI is used. The radial cuts marked on the 
image represent the approximate slit positions of the observations. 
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they are : weakening of the polarization signal due to scattered light within the in strument 
and fr om nearby quiet Sun; peculiarity of the atmosphere such as discussed by iHanaoka 
( I2OO5I ): the methods used to derive LOS magnetic field (inversion of Stokes profiles is 
yet to be established); inclination of t he magnetic field; and most iniporta ntly the height of 
formation which is highly ambiguous (jSocas-Navarro fc Uitenbroekll2004l ). We will discuss 
these issues in detail after looking at the results from the analysis of velocity gradients both 
at the photosphere and chromosphere. 



5.2. Velocity Gradients and the Nature of the Magnetic Fields 

The velocities calculated using COG method show a typical behavior of Evershed flows 
both at the photosphere and chromosphere. That means the limb side penumbra shows 
red shift with respect to quiet Sun where as center side penumbra shows blue shift at the 
photospheric heights. The situation is exactly opposite at chromospheric heights which are 
consistent with the well known inverse Evershed effect. COG velocities in umbral regions 
both at photosphere and chromosphere are smaller compared to penumbral regions. 

The plots of bisector velocity differences (Eq. [5]) as a function of photospheric magnetic 
field strength are shown in Fig. H] for the Fe I line. Top panel in this figure is for all the 
points over the total field of view (FOV), where as the bottom left and right panel is for 
umbral and penumbral regions of the observed spot, respectively. Note from these figures 
that the large number of points correspond to umbra have smaller velocity gradients than 
the penumbra. Closer investigation of bisector velocity differences of Fe I line show the flow 
pattern consistent with the well known Evershed flow. That means larger portion of the 
limb side penumbra shows net downfiow (both core and wing side bisectors show redshifts) 
and disk center side penumbra shows net upfiow (both core and wing bisectors show blue 
shifts with respect to the refe rence). For the definition of the net up- and down-flow see 



Balasubramaniam et al.l (119971 ) 



We also found that the bisector velocity gradients and COG velocities observed in Fe I 
show a good correlation in agreement with the earlier observations. That means, larger the 
COG velocity larger the velocity gradient at the photosphere. 

In Fig. Othe plots of bisector velocity gradients for v/s the photospheric magnetic 
fields are shown. Top panel in this figure is for the total FOV and left bottom panel for 
umbral region and right bottom panel for penumbral region. Notice the increase in bisector 
velocity gradients measured in Ha with increase in the photospheric magnetic field strength. 
Plots of velocity gradients along two radial slices of sunspot (marked as 1 and 2 in Fig. 
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ls]) also show that the larger velocity gradients at the chromosphere are located at which 
strong photospheric fields are observed (Fig. [6]) . Top panels in Fig. [6] show the plots of 
photospheric field strengths along the radial cut 1 and the radial cut 2. The corresponding 
plots of velocity gradients are shown in the bottom panels. Except at few locations (mostly 
in the limb side penumbra such as shown in the bottom right panel of Fig. [6]) which show 
less variation in the values of velocity gradients inspite of variation in the photospheric field 
strengths (which is evident also in Fig. [5]), most of the places the velocity gradients are 
larger where the photospheric magnetic field strengths are larger. 

Closer examination of the bisector wavelength positions of Ha Stokes / profiles in dif- 
ferent regions of the sunspot with respect to the reference wavelength reveals the following 
results (see figure [7]). 

• In the limb side penumbra both line core and wing side bisectors show blue shift with 
respect to quiet Sun. The shifts in the fine core bisectors are large compared to the 
line wing bisectors indicating the net upfiow. 

• In the umbral region also both line core side and wing side bisectors show blue shifts 
with respect to the reference wavelength. The shifts in the line core side bisectors for 
the umbra are almost comparable to that of the limb side penumbra. However, the 
shifts in the line wing side bisectors for the umbra are much smaller with respect to 
that of the limb side penumbra. This would indicate again the net upfiows but with 
larger velocity gradients. 

• In the center side penumbra the line wing bisectors show redshift in most of the places 
with respect to the quiet Sun. At few locations they show a small blue shift or no shift. 
The line core side bisectors show slight blue shift with respect to quiet Sun reference 
wavelength position. 

It is found through these analyses that the velocity gradients are larger in the umbra 
at the chromospheric heights compared to the penumbra (see Fig. [5]). Which is exactly in 
contrast with the fiow pattern observed at the photosphere. The analysis of bisector velocity 
differences also indicate accelerated upfiows in the umbral region. 

There is a little concern due to the residuals present after the Co I line blend is removed. 
To make sure that there are no artifacts introduced due to the residuals, the bisector velocity 
differences calculated by considering the spectral region of the line which is not affected by 
the blend are shown in Fig. [8] as a function of photospheric field strength. This figure also 
shows the trend that the velocity gradients increase with increase in photospheric magnetic 
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field strength confirming the observations made with the blend removed full intensity profiles. 
As expected, the gradients are smaller due to the smaller wavelength regions considered in 
this case. 

To summarise, wavelength shifts analysis of Stokes I bisectors show that the velocity 
gradients are larger in the umbral region than in the penumbral region at the chromospheric 
heights. Most importantly, accelerated upfiows are observed in the umbral region. In other 
words, LOS velocity increases upward more rapidly in the umbral region than in the penum- 
bral region at the chromosphere. 

Let us now address some of the possibilities, mentioned in the beginning of this section, 
which can cause the observed weaker chromospheric field in the umbral region from the point 
of view of the velocity gradients. 

We found that the stray light within the instrumen t is less than 2% by comparing the 



quiet Sun spectrum with the atlas (IWallace et al.l 120001 ). This amount of stray light is too 



small to account for the observed weaker Stokes V signals in the umbral region. Also, it can 
not account for the observed velocity gradients in H^. 

The next question raised was about the reliability of the methods used to estimate the 
LOS magnetic field. Despite the availability of more accurate estimation of photospheric 
magnetic field through inversion, we have used the COG method to maintain the uniformity 
while comparing the LOS field at the photosphere and chromosphere (Fig. [2]). This is 
because the inversion of Ha profiles to derive magnetic field is yet to be estabhshed. As 
discussed in section WA\ the underestimation of the LOS field obtained from Fe I is expected 
to be more while the values obtained from Ha will be less as confirmed through numerical 
simulations. Also, the observed umbral fields are larger at the photosphere where as they are 
smaller at the chromosphere. Hence the departure of LOS magnetic fields from the actual 
values due to the method used to estimate them can not explain the observed weakening of 
the magnetic fields. 

One of the possibilities which can cause the r eduction in polarization (in effect the 



weaker magnetic field) proposed by iHanaokal (120051 ) is the peculiarity of the atmosphere. 



However, the kind of peculiarity discussed by him can not explain the velocity gradients and 
the Stokes profiles of Ha discussed in this paper. 

It is a general wisdom gained from the extrapolation technique that the magnetic topol- 
ogy at the chromosphere is very much dependent on the field configuration at the photo- 
sphere. Extrapolation techniques like potential field approximation suggest that larger the 
field inclination at the photosphere chromospheic fields should also show larger inclination. 
From the inversion of Fe I line it was found that the fields are oriented more close to LOS 
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in the umbral region than in the penumbral region. Hence we expect larger field orientation 
in penumbral region at the chromospheric heights also. As mentioned in section 14.21 
shows Stokes Q and U signals which are above the noise level in penumbral region but not in 
umbral region indicating that the orientation of the fields are larger in the penumbral region 
compared to the umbral region even at the chromosphere. Hence we believe that the weaker 
LOS field observed in the umbral region may not be due to the larger orientation angle with 
respect to the LOS. This scenario is also verified using the velocity gradients observed at 
the chromosphere, as larger inclination means smaller LOS velocity in contradiction to the 
observed velocities (Fig. [5]). 

Another major difficulty in interpreting the obse rvations is the ambiguity in its height 



of formation. Studies based on the response functions by lSocas-Navarro fc Uitenbroekl (120041 ) 
show that the Ha is sensitive mostly to chromospheric magnetic field in the umbral model. 
While in the quiet Sun model it shows sensitivity to both photospheric and chromospheric 
magnetic fields. Since the magnetic fields at the photosphere is large, the magnetic field 
measured in quiet Sun regions will mostly be photopsheric. If the Ha were to show the 
sensitivity to magnetic field in the penumbral model similar to that of quiet Sun model then 
one would expect larger field values in the penumbra which will be an average of photospheric 
and chromospheric fields while the umbral fields are exclusively chromospheric. To confirm 
this, more forward modeling is needed which includes the comparison of Ha line formation 
at different regions with different field configurations. More over, response of Ha line to 
various physical parameters needs to be studied and consistently explain the observations 
such as velocity and velocity gradients presented in this paper. 

Other possibility which can cause the observed weaker chromosp heric fields in the umbral 



region is the faster divergence of stronger fields as suggested by iBalasubramaniam et al. 



( I2OO4J ). This scenario may consistently explain both the observed magnetic properties (Fig. 
El) as well as the velocity properties (Figs. [S)/1H])- Because of the decrease in the field strength, 
the plasma becomes less inhibited by the magnetic fields and hence more free to move. That 
means if the stronger fields diverge faster compared to weaker fields then the the velocities 
should increase faster in the stronger field regions which are consistent with the observations. 



5.3. Stokes V Amplitude Asymmetry 

From the close examination of Stokes V profiles of Ha we found that the Co I line has 
not intruded up to the extent that the amplitudes are affected. Hence the results presented 
from Ha based on amplitude asymmetry are reliable. 
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From the Fig. [9] it is clear that the amphtude asymmetry of Fe I hne tends toward 
zero with increase in photospheric magnetic field strength. This means, at the photospheric 
heights the gradients are smaller in umbra which corresponds to the region of strong fields. 
In contrast amplitude asymmetry observed in tend to increase with photospheric field 
strength as shown in the Fig. [TOl This implies that the field gradients are larger in the umbra 
at the chromosphere compared to penumbra. Hence the analysis of amplitude asymmetry 
confirm the results obtained from the bisector analsysis in section 15.21 



6. CONCLUSIONS 

In this paper we have used the multiwavelength spectropolarimetric tool to understand 
the stratification of the magnetic and velocity fields in the solar atmosphere. Out of the two 
lines considered for spectropolarimetry, one forms at photospheric height (Fe I A6569) and the 
other spans almost from the upper photosphere to upper chromosphere {Ha). Hence, these 
lines were useful in studying the connection between physical parameters at the photospheric 
and chromospheric heights. The main physical parameters studied in this paper are the 
magnetic and velocity fields in an active region. 

As discussed in section 15. Ij the LOS magnetic field measured in chromosphere is in 
general weak compared to its photospheric counterpart. The weakening of the chromospheric 
field is much faster for the corresponding strong photospheric field. The magnetic field 
strengths observed in the umbral chromosphere are much weaker than those expected from 
the extrapolation of the photospheric magnetic field. For instance, the field strength inferred 
through Ha observations is about 400 G where as the field strength obtained through the 
extrapolation of the observed photosperic field under potential field approximation to an 
height of 2000 km is about 1000 G (assuming that the height of formation of Ha is about 2000 
km). Various possibilities have been discussed which can cause the weaker fields observed in 
the umbral chromosphere. The most probable ones are the fast divergence of the stro nger 



fields when they propagate upward in the atmosphere (IBalasubramaniam et al.ll2004l ) and 



ambiguity in the height resolution of Ha magnetic sensitivity which may be photospheric 



and / or chromospheric depending on the region of observation (jSocas-Navarro fc Uitenbroek 



20041 ). If former is the reason, then it can explain the observed properties of both velocity 
and the magnetic fields presented in this paper. This is because, as shown in Fig. [31 observed 
chromospheric fields are systematically weaker at the locations where strong photospheric 
fields are observed. If the weaker field strengths observed at the umbral chromosphere are 
truly of solar origin then this implies that the umbral fields decrease more rapidly with 
height compared to penumbral fields. Rapid decrease in field strengths cause rapid increase 
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in velocity with height as there is a small upflow in the umbral photosphere. Observations also 
show that the velocity increases more rapid ly in the umbral region cona pared to penumbral 
region (section [5.21) . Earlier observations by lGosain &: Choudharyl (120031 ) have also indicated 
the weakening of the magnetic field which is larger for stronger fields. In their observations 
quicker weakening of the stronger fields is not apparent, probably, because of the lines (Mg I 
bl and b2 at A5173 and A5184) used to infer chromospheric magnetic field that originate at 
the lower chr omosphere. While the ob s ervat io ns based on if presented in this paper as well 
as earlier by iBalasubramaniam et al.l (120041 ) ; iHanaokal (120051 ) show clearly that the strong 
fields weaken quickly with height because in these works, to infer chromospheric magnetic 
field only the spectral region close t o its liiie core is considered which samples mostly the 
higher layers of the chromosphere (iRuttenI 120071 ). Hence there is a possibility that the 
weaker LOS field strengths observed in the umbral chromosphere are caused due to the 
faster divergence of the stronger fields. However, we would like to caution that the reliability 
of the COG method in estimating the LOS magnetic field strengths discussed in this paper 
(section 14.11) is for simple solar atmospheric model. But, in reality the solar atmosphere 
may be complicated. More studies on the magnetic and velocity response functions for 
Ha in different regions will help in better interpretation of the observations. Simultaneous 
multiline spectropolarimetry, which includes Ha and at least one more line which is formed 
at the chromosphere (preferably Infrared Ca Triplet line at A8542) and a photospheric line, is 
needed to get further insight into the physical processes that take place in the chromosphere. 
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Fig. 1. — Typical Stokes / profile with the bisectors marked in diamond symbols. Dotted 
line in the top panel shows the spectral region of the profile after removing the blend. 
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Fig. 2. — Scatter plot of photospheric and chromospheric LOS magnetic field derived using 
COG method applied to Fe I(A6569) and respectively. Solid line is a linear fit to the 
field strengths correspond to the penumbral region. 



Fig. 3. — On the right side of this figure shown are the plots of LOS magnetic field strengths 
at the photosphere (dashed lines) and the chromosphere (dotted lines) along two radial 
slices of the sunspot. For reference, these two radial cuts are marked on the SOHO/MDI 
intensitygram. The arrow on the sunspot image indicates the disk center direction. 
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Fig. 4. — Plots of velocity gradients at the photosphere v/s photospheric magnetic field 
strength. 
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Fig. 5. — Plots of velocity gradients at chromosphere derived using Ha with the bisectors 
considered for the full line profile v/s the total field strength (photospheric) . Solid lines are 
the linear fit to the data points. 
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Fig. 6. — The top panels in this figure show the plots of magnetic field strength at the 
photosphere along two radial slices of the sunspot. The radial slices are same as those shown 
in Fig. [3l The corresponding plots of velocity gradients at the chromosphere are shown in 
the bottom panels. 
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Fig. 7. — Representative bisectors of Ha line profiles correspond to different regions of the 
sunspot are plotted on an average quiet Sun profile. The solid line represents the reference 
wavelength (which is a COG wavelength position of the nearby average quiet Sun profile). 
Dotted line represents the bisectors location typical of umbral profiles, dashed line-limb side 
penumbra, dash dotted line- center side penumbra. 
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Fig. 8. — Same as Fig. [5] but with the hmited spectral range about the hue core considered 
for calculating velocity gradients to avoid the influence of blending due to Co I line. The axes 
scales are kept same as Fig. [5] to indicate that the velocity gradients are smaller compared 
to the case when the full Stokes / profile is considered. 
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Fig. 9. — Plots of amplitude asymmetry of Stokes V profiles of Fe I line v/s the total field 
strength (photospheric) . The sohd lines are the hnear fits to the data points. 



-27- 



full FOV 



-0.2 



500 1000 1500 2000 2500 

B(Photosphere) Gauss 



umbrol region penumbrol region 




1000 1500 2000 2500 500 1000 1500 

B(Photosphere) Gauss B(Photosphere) Gauss 



Fig. 10. — Plots of amplitude asymmetry of Stokes V profiles of line v/s the total field 
strength (photospheric) . The sohd lines are the hnear fits to the data points. 



